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The caloric curve for mononuclear configurations is studied with a schematic model. We inves- 
tigate the dependence of the entropy on the density and effective mass profiles. A plateau in the 
caloric curve is a direct result of decreasing density and the destruction of correlations with in- 
creasing excitation. The mononuclear regime is metastable with respect to binary fission at low 
excitation energy and unstable with respect to multifragmentation at high excitation. The statisti- 
cal framework presented here is suitable to treat scenarios where experimental conditions are set to 
favor population of highly excited mononuclear states. 
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Heavy compound nuclei (CN) are metastable objects; 
local mononuclear entropy maxima separated from the 
more stable dinuclear states by transition states of signif- 
icantly lower entropy. Standard statistical-model treat- 
ments of CN decay are predicated on a time-scale sep- 
aration between the CN formation time, and the time 
scales for simple (mostly single-particle) decay modes 
and the massively collective decay processes we call fis- 
sion. The former are usually treated using the prescrip- 
tions offered by WeisskopfyLj or Hauser and FeshbachQ 
while the latter are usually treated by a transition-state 
formalism initially developed for chemical reactions by 
EyringQ. The distinction between decay modes can be 
bridged in concept^] and practice^ however the mas- 
sively collective decay channels can be retarded by tran- 
sient delays 0- 

It has been known for almost 4 decades that below ap- 
proximately 1/3 of the saturation density p Q , a- matter 
has a lower free energy than uniform nuclear matter • 
However mononuclear configurations at reduced density 
can be reasonable subjects for statistical decay treat- 
ments as long as they are either metastable or protected 
from massively collective decays modes by transient de- 
lays. Specifically with increasing excitation energy, an 
equilibrium (i.e. local maximal entropy) mononuclear 
density profile can be reached, the decay of which can be 
treated with only minor modifications to the well known 
formalisms, as long as a time-scale separation exists be- 
tween formation, equilibration and all conceivable decay 
modes. 

This work does not deal with the important issue of 
time-scale separation 9] , only with the gross statistical 
properties of mononuclei at high excitation energy. In 
particular, we show that the relaxation of the density pro- 
file of mononuclei, in pursuit of maximal entropy, causes 
the Caloric curve, T{e) (where T is the statistical nu- 
clear temperature and e is the excitation per nucleon) 
to flatten out and exhibit a quasi-plateau. We believe 
this explains the nature of the Caloric curve first studied 
by Wada et al.[lO( and later by Pochodzalla et al.[ll|, 
and for which systematics have recently been analyzed 
in detail by Natowitz et al. |T^ . 



The approximate saturation of the statistical temper- 
ature is primarily due to density reduction, but is also 
influenced by the evolution of the effective mass of nu- 
cleons in the nuclear medium. The first effect is just the 
sequestration of energy in the potential energy of nuclear 
expansion. The energy spent on expansion reduces the 
thermal part of the total excitation in much the same 
way as the collective rotational energy does in the case 
of high angular momentum CN. 

The ratio of the effective mass to the bare nucleon 
mass m* /m differs from 1 due to the finite range of the 
nuclear force and the time non-locality of the interac- 
tion. The former effect, which is responsible for making 
the optical model potential energy dependent, reduces 
m* jra by a density dependent factor m%{p) which must 
return to 1 at low density. The time nonlocality can be 
thought of as the coupling of low -lying s urface modes to 
single-particle degrees of freedom [l3| [ li j . This collective 
effect brings strength down from high energy, increasing 
the many-body density of states at low excitation energy. 
The effective mass factor, m w (p',T), accounting for this 
relocation of strength, while greater than 1 at low energy 
and localized on the surface of the quantum drop, must 
return to 1 in the limit of high excitation as well as low 
density gradient. 

We confine our analysis to a one parameter descrip- 
tion of expansion and literature descriptions of how 
the effective-mass terms evolve with density and exci- 
tation. Our approach combines the physically transpar- 
ent picture of maximal-entropy mononuclear configura- 
tions found in the recent work by Toke et al. |15j . with 
effective-mass logic similar to that found in the works of 
Natowitz, Shlomo and collaborators |l6j . 

A quantum drop of degenerate Fermi liquid has an 
entropy given bv[l7|. 



S M = 2VuU = 2^a{E* T - E c ) = 2^/aA(e - e c ), (1) 

where a is the level-density parameter and the thermal, 
total and compressional energies are U, Ej,, and Ec, re- 
spectively. With total particle number A, e and ec are 
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the total and compressional energies per particle. In the 
local density approximation f LDA) [lq] . the level density 
depends on the nuclear profile, the local Fermi momen- 
tum k and the effective mass 191 201, 



Pr( r ) 



[h^fc2(r)/2m* 



-dr. 



(2) 



The density profiles p T (r) of the two isospin partners 
(with index r) are taken to be the same functional form, 
scaled in proportion to the number of nucleons. The 
native (e = MeV) radial profiles are of the "standard" 
type with a Gaussian derivative, 



P„(r) = §(l 



erf (- 
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with effective sharp radius R — TqA 1 ^ (r a = 1.16 
fm) and surface width b = 1.0 fm. The expansion is 
limited to the one-dimensional self-similar family, i.e. 
p(r, c) = c 3 p n (c* r). The expansion parameter c is found 
by maximizing the mononuclear entropy, 



dS M (E* T ,p(c)) 
dc 1 



0. 



The collective energy involved in expansion is taken 
as the simple parabolic form, involving only the central 
density, suggested by Friedman^, E c — /l*e;,(l — ^p-) 2 ■ 
We have used Eb = 8 and 6 MeV (red a nd g reen curves 
in figures 2-4) in the present calculations 22J. 

Execution of eq. 4, not only finds the metastable 
mononuclear expansion but also insures that the surface 
pressure is zero. This procedure is therefore logically dif- 
ferent from the physically unreal, but true equilibrium 
condition found by placing a drop in a box and having a 
surrounding vapor supply a pres sure. 

Following Prakash et al. j2(j, and De et al. [2^| we 
choose a phenomenological form for the effective mass 
that captures the underlying science, 



!* =imu )[ mu ] = {l-aM )[1 _ m m ] 
with 



(5) 



a = 0.3, (3{T) = 0.4A 1 / 3 exp[-(TA 1 / 3 /21) 2 ]. (6) 

(The T dependence in (3{T) requires performing an it- 
eration of the outlined procedure. This modification is 
straightforward and will not be described.) The effective- 
mass factor is suppressed in the bulk, peaks at the 
surface [2{| and degrades to 1 with decreasing density and 
increasing thermal energy. These two many-body effects, 
to a large extent, offset one another in near ground-state 



nuclei, yielding a ~ ^4/8.6 for unexpanded 197 Au, the 
nucleus considered here. However the destruction of the 
cooperativity encoded in these two effective-mass terms 
does not occur on identical energy scales. While the de- 
tailed density and excitation energy dependence of these 
terms is unknown, the present work shows how the gross 
effects captured by these terms couple with expansion to 
dictate the form of the Caloric curve. 

The excess entropies, S M — Sm — Sm(c — 1), are 
shown in Fig. 1. The maximum entropy determines the 
equilibrium expansion parameter (c) as a function of e. 
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(4) Fig. 1 : Representative calculation of excess entropy as a 

?m( 



function of expansion: Sf,j(c) = Sm ~ Sm(c = 1). 
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The reduction of central density with excitation p c (s) 
is shown in Fig. 2 without [to* = 1] (black) and with 
[m*/m = mk{p)m ul {p' ,T)\ (red: = 8 MeV and green: 
Eb = 6 MeV) effective mass considerations. (Consider- 
ation of TOfc (/?) alone exhibits a reduction in the central 
density similar to that with •mk{p)m UJ {p' ,T), while con- 
sideration of mk{p)m UJ {p r ) leads to approximately the 
same p c (e) as with to* /to = 1.) Without the effective- 
mass terms, the p c (e) dependence is almost identical with 
the (extended) finite-temperature Hartree-Fock calcula- 
tion reported in [l^. As one should expect, the de- 
crease in density is more substantial with the smaller 
compressibility. The central densities implied from the 
maximal-entropy procedure used here with to* /to = 
mk(p)m LJ (p' ,T) and Sb — 6 MeV are similar to those 
extracted from caloric curve data in [l6| (diamonds) . On 
the other hand, p c (s) does not drop quickly, over a nar- 
row a range of energies, as is suggested by the analysis 
of Coulomb barriers by Bracken et al. [2]} (circles). 

Figure 3 displays Sm for the native density profile (dot- 
ted), as well as Sm and S M for to* /to = 1 (black), 
to* /to = m k (p)m UI (p' ,T) for e b = 8 MeV (red), and 
£6 = 6 MeV (green). The excess entropy allowed by re- 
laxing the density profile (black dashed) is evident as is 
a reduced rate of entropy growth at for e < 4 MeV when 
the effective mass is modelled by equations 5 and 6. This 
"reduction", is due to decreasing to* /to in the surface 
region with increasing e. 
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Fig. 2: Central densities for 197 Au with [m* = 1, e b = 8 

MeV] (black), [to* = m k {p)m^{p' \T), e b = 8 MeV] 
(red), and [to* = mk{p)m iA] (p l \T) , e b = 6 MeV] (green). 
Also shown are the densities extracted from apparent 
level density parameters (diamonds) 1261 and from 
Coulomb barriers (circles) |27fl. 
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Fig. 3: Sm for a native density profile (dotted), and Sm 
(solid) and S M (dashed) for the cases shown in Fig. 2. 

The statistical temperatures, found by differentiation, 



T = A/[—) v , (7) 

are shown in Fig. 4. (The nuclear volume is not a ther- 
modynamic control variable, rather a volume sufficiently 
large to make the pressure zero is. This volume is un- 
affected by the changes in the nuclear volume. From 
another perspective, as the external pressure is zero, no 
work is done by the nuclear expansion.) As shown in |l5j 
and inferred in |16[ . the relaxation of the density profile, 
substantially flattens the temperature rise with e (com- 
pare black dashed and solid curves) . The inclusion of the 
full effective mass form (function of p, p', and T) increases 



T for e < 4 MeV and decreases it for e > 5 MeV (com- 
pare solid black to colored curves), changes that give the 
appearance of a plateau. Decreasing the compressional 
energy constant (red to green), reduces the value of the 
temperature of the pseudo-plateau. 
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Fig. 4: T for the same cases as shown in the previous 
two figures. 

Thermodynamic stability of a single phase requires 
the total St{E) have a positive first derivative (positive 
temperature) as well as being concave (implying a posi- 
tive heat capacity C v as (§]^)v = —l/T 2 C v ). However 
the same cannot be said of the constrained mononuclear 
Sm(E). Comparison of Sm{E) to Smf, the entropies for 
multifragmentation, is complicated by the prescription 
for the free volume needed for the latter [2^] [2^]. Never- 
theless one can compare SM(m*/m = 1), at extracted 
temperatures, to Smf &s a function of volume. Taking 
the latter from the ideal phase space model of Das Gupta 
and collaborators 30] indicates: 1) The volume capturing 
99.75% of the mononuclear matter is much smaller than 
reasonable volumes for multifragmentation. 2) Smf will 
exceed Sm somewhere above s > 3 MeV if a freezeout 
volume of 3 times the unexpanded volume is used. This 
crossing implies that the mononucleus becomes unstable 
with respect to the multifragmented state trapped in a 
hypothetical freezeout volume. As the latter occupies 
a much larger volume, the multiparticle phase space re- 
gions are divorced indicating that transient times, as well 
as pure statistical rates, are needed in a decay treatment. 

An irreversible phase change (finite free energy change) 
from an initially populated metastable or kinetically 
trapped mononucleus to the higher entropy multifrag- 
mented state would proceed with evolution of the kinetic 
energy and thus event averages, sampling a time distri- 
bution, would yield fluctuations in excess of those ex- 
pected for an equilibrium system. We therefore sug gest 
that excessive fluctuations of the kinetic energy |3j], re- 
sult from an irreversible transition from a lower entropy 
(larger free energy in a canonical treatment) but kineti- 
cally trapped, mononuclear phase to the multifragmented 
phase. This is equivalent to arguing that the ensemble 
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(the event sample) is non-ergodic. However this tran- 
sition is not required to explain a plateau in a caloric 
curve, nor are large fragment multiplicities or significant 
final-state correlations. 

As concerns fragment multiplicity, the characteristic 
(evaporative) decay time r ev — £, with the decay width 
calculated using realistic, but simple forms for the neu- 
tron decay width F ra [33, is of the order of 10-15 fm/c 
across the plateau. In a standard (binary transition- 
state) statistical model, many fragments would be evap- 
orated as the CN retraces the Caloric curve going to the 
left, contracting as the sequence continues. As the tem- 
perature does not change significantly across the plateau, 
the issue of sequence is irrelevant. 

The relevance of such a statistical model depends on 
there being a time-scale separation between mononu- 
cleus formation r f or and the decay times associated with 
mononuclear r ev and multifragment T m f final states. It is 
unlikely that T/ or < r ev in heavy-ion reactions, but this 
condition might be met for light ion, it and p induced re- 
actions. The time scale relevant for T m t is the maximum 
of a purely statistical time, related to the entropy of the 
transition state Smf, and the delay time associated with 
the amplification of density fluctuations. Work on the 
difficult issue of the spinodal decomposition time scale0 
suggests that r ev < T m f in the plateau. A clear analogy 
between the transient delay for fission 6] and the ampli- 
fication time for density fluctuations leading to spinodal 
decomposition can be made. 

The present work indicates that a purely statistical 



treatment of the decay of a highly excited nucleus (pre- 
suming formation of an equilibrated mononucleus) would 
be one which evaporates from metastable (local maxi- 
mal entropy) expanded systems and has a decay branch 
with transient delay for spinodal decomposition. The 
significance of the delay time is that even if Smf > Sm 
the decay might still be determined by the mononuclear 
phase space. We note that the logic for treatment of 
the evaporative-decay channels outlined above is sub- 
stantially different from that offered in the "Expanding 
Emitting Source" model proposed by Friedman[2l|. In 
EES, the expansion is dynamically produced rather than 
the result of maximizing the entropy. 

This work makes use of several simplistic assump- 
tions. The LDA is known to yield inaccurate level density 
results 18] and the self-similar expansion is a restriction 
that limits the entropy. The expression for the entropy 
is also simplistic in that it ignores the continuum and 
all detailed quantum structure. Nevertheless, it seems 
that a near plateau of the Caloric curve is an unavoid- 
able consequence of an expanding metastable Fermion 
system (with a finite-range force) and the resultant evo- 
lution (destruction) of collectivity with expansion and 
excitation. 
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